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Abstract 

We review recent progress at the two e'’'e“ B factories. The first 
measurement of CP violation and the prospects for measuring all the 
angles of the unitarity triangle are discussed. 


1 High Luminosity e+e B factories 

1.1 Introduction 

Many important discoveries in high energy physics were announced in the 
Rencontres de Moriond series. For example, at the 1987 Rencontres de 
Moriond in a presentation by H. Schroeder, the ARGUS collaboration de¬ 
scribed the surprising discovery of large Bd — Bd mixing|]I|. The implications 
of this observation were important for the experimental program on CP viola¬ 
tion. As a result, the theoretical proposal of Bigi, Carter and SandaQ] that 
there would be large CP asymmetries in hadronic B decays to CP eigen¬ 
states became experimentally feasible. This initiated the program which 
would eventually lead to the plans for the high luminosity e’''e“ B factories. 

Two international high luminosity B factories have now been com¬ 
pleted and are collecting data in the United States (PEPII/BABAR) and 
Japan (KEKB/BELLE) P],^. The goal of these experiments is to measure 
CP asymmetries in the B system and test the foundations of the Standard 
Model. 
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The accelerators were commissioned with remarkable speed starting in 
late 1998. The experiments starting physics data taking in 1999. By the 
time of the XXXVI Rencontres de Moriond in winter 2001, there were hints 
(at the 1.7 ct level) of CP violation from both BABAR and BELLE. In the 
summer of 2001 , these two experiments announced the observation of the first 
statistically significant signals for CP violation outside of the kaon system. 

1.2 CP Violation Induced by Mixing 

In 1973, Kobayashi and Maskawa (KM) first proposed a model where CP vio¬ 
lation is incorporated as an irreducible complex phase in the weak-interaction 
quark mixing matrix |^. The idea, which was presented at a time when only 
the M, d and s quarks were known to exist, was remarkable because it re¬ 
quired the existence of six quarks. The subsequent discoveries of the c, b and 
t quarks, and the compatibility of the model with the CP violation observed 
in the neutral K meson system led to the incorporation of the KM mech¬ 
anism into the Standard Model, even though it had not been conclusively 
tested experimentally. 

The unitarity of the CKM matrix implies that the existence of three 
measurable phases. In the “Nihongo” convention, these are denoted 



while at SLAC these angles are usually referred to as /3, a and 7 , respectively. 
A non-zero value of (pcp results in the time dependent asymmetry 


, R{B^ ^ fcp) - R{B^ ^ fcp) 

^ ^ fcp) + R{B^ ^ fcp) 


^fsm2(j)cp ■ sm{Am ■ (t 2 ±ti)), ( 2 ) 


where is the CP eigenvalue (±1), Am denotes the mass difference between 
the two mass eigenstates and ti and ^2 are the proper time for the tagged- 
B and CP eigenstate decays, respectively. The -|- sign corresponds to the 
case where the B^ and B^ are in an even L orbital angular momentum state, 
the — sign obtains for odd L states. A determination of Af thus provides a 
measurement of sin 20 cp. 

We note that due to the restrictions of quantum mechanics, time inte¬ 
grated asymmetries at the T (45*) resonance (which corresponds to the — sign 
in the above equation) are identically zero. Therefore, one must make time 
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dependent measurements. Since the pairs of B mesons are produced nearly 
at rest in the usual arrangement at threshold, the T (dS*) center of mass frame 
must be boosted. This is accomplished by the use of beams with asymmetric 
energies. For example at KEK-B, /dy ~ 0.43, (at PEP-II, /dy ~ 0.56) and the 
typical B meson decay length is dilated from ~ 20/im to ~ 200/im, which 
is measurable with solid-state vertex detectors close to the interaction point. 
At the symmetric CESR storage ring which uses equal energy beams, these 
measurements are not possible. 

1.3 Accelerators and Detectors 

The e+e“ B factory experiments require substantial improvements in accel¬ 
erator technology to achieve luminosities in excess of 3 x 10^^/cm^/sec and 
data samples with integrated luminosities of 0(100 fb~^) as well as some 
improvements in detector technology. The properties of the B factory accel¬ 
erators are summarized in Table 1. The properties of the individual bunches 
are not dramatically different from what had been achieved at conventional 
storage rings such as CESR. Instead to obtain high luminosity, the asymmet¬ 
ric B factories use double storage rings and very large numbers of bunches. 
This is a domain of accelerator operation which had not been fully explored 
previously and to avoid instabilities involving coupling of bunches requires 
the operation of complex feedback systems. 

To minimize parastic collisions between incoming and outgoing bunches, 
KEK-B chose the crossing angle approach while SLAG takes advantage of 
the energy asymmetry and uses magnetic separation to achieve this sepa¬ 
ration. A small crossing angle of ±2.1 mrad had been successfully used at 
CESR to achieve high luminosity. KEK-B employs a somewhat larger ±11 
mrad crossing angle. There is still some worry that coupled oscillations of 
the longitudinal and transverse degrees of the beam (synchro-betatron os¬ 
cillations) may be excited, so KEK-B continues R±D to develop special RE 
cavities to rotate the bunches so that they collide head on at the interaction 
point. The magnetic separation approach gives rise to somewhat larger ex¬ 
perimental backgrounds and imposes tight constraints on the engineering of 
the interaction region. 

The performance achieved by both accelerator complexes in a very short 
time is quite remarkable. Both started physics running in 1999. By the end of 
the summer of 2001, PEP-II had integrated about 40 fb“^ while KEK-B had 
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integrated 33 fb“^. The peak luminosity for PEP-II was 3.9 x 10^^/cm^/sec 
while KEK-B had achieved 4.5 x 10^^/cm^/sec. In the case of KEK-B, the 
beam currents are still far below the design values. (Also note that the 
operating parameters currently used are constantly changing and somewhat 
different from those in the design (Table 1)). 

The experimental measurements of CP violation require tagging of the B 
flavor, measurement of the time dependence of the decay and high momentum 
particle identihcation. As shown in Table 2, in contrast to the accelerators, 
the detectors are quite similar with the exception of the particle identihcation 
systems §,0. 


Table 1: 

Design Parameters for B-Factory Accelerators 



KEK-B 

SLAC 



LER HER 

LER 

HER 

Energy 

E(GeV) 

3.5 8.0 

3.1 

9.0 

Luminosity 

L(cm“^s“^) 

1 X 10^^ 

3 X 10 

33 

Collision mode 


±llmr (crab) 

Head-on 

Circumference 

C(m) 

3018 

2199 


Beta Function 


100/1 100/1 

37.5/1.5 

75/3 

Tune shift 


0.05/0.05 

0.03/0. 

03 

Emittance 

e^/ey{n-m) 

19/0.19 19/0.19 

64/2.6 48.2/1.9 

Energy spread 

aE/E{lQ-^) 

7.7 7.2 

9.5 

6.1 

Total Current 

J(A) 

2.6 1.1 

2.14 

0.98 

No. of Bunches 

Nb 

5120 

1658 


Bunch spacing 

S'ij(m) 

0.6 

1.26 


RE frequency 

/rf(MHz) 

508 

476 


RE voltage 

Rc(MV) 

22 48 

9.5 

18.5 

Cavity type 


ARES super 

1 -cell normal 

No. of Cavities 

Nc 

28 60 

10 

20 


The most challenging requirement for the detector is the detection of short 
lived vertices. The experiments at both B factories use double-sided silicon 
strips to deal with this requirement. These have been used previously at 
ALEPH, DELPHI, and at CLEO II and III. Nevertheless, the implementation 
is still far from routine. 

The other difficult requirement for the detector is the separation of kaons 
from pions. At high momentum, this is needed to distinguish B^ —> tt+tt” 
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from K~n^. The former is a spectator b ^ u mode which is used to 

determine one of the angles in the unitarity triangle while the latter is an 
example of a gluonic penguin b sg* decay. Similarly, the experiments must 
distinguish the Cabibbo suppressed mode B~ D^K~ from the much more 
abundant and kinematically similar Cabibbo favored mode B~ D^t[~. At 
lower momenta, particle identihcation is essential for flavor tagging. 

Two approaches for high momentum particle identihcation have been 
implemented. Both are based on the use of Cerenkov radiation. At BELLE, 
aerogel Cerenkov radiators are used. Blocks of aerogel are readout directly 
by hne-mesh phototubes which have high-gain and operate comfortably in 
a 1.5 Tesla magnetic held. Since the threshold for the aerogel is around 
1.5 GeV, below this momentum K/tt separation is carried out using high 
precision time-of-hight scintillators with resolution of 95 ps. The aerogel 
and TOE counter system are complemented by dE/dx measurements in the 
central drift chamber. The dE/dx system provides K/tt separation around 2.5 
GeV in the relativistic rise region and below 0.7 GeV. For high momentum 
kaons, an efficiency of 88% with a misidentihcation probability of 9% has 
been achieved. 

At BABAR, Cerenkov light is produced in quartz bars and then trans¬ 
mitted by internal rehection outside of the detector through a water tank 
to a large array of phototubes where the ring is imaged. The detector is 
referred to by the acronym DIRC. It provides particle identihcation over the 
full momentum range for particles that are energetic enough to reach it. Ad¬ 
ditional particle identihcation is provided by dE/dx measurements from the 
drift chamber and from the 5-layer silicon detector system. 

For tracking, both B-factory experiments use conventional drift cham¬ 
bers with a low-Z helium based gas to minimize multiple scattering and SR 
(synchrotron radiation) backgrounds. 

To detect photons and electrons, the B factory detectors use large arrays 
of CsI(Tl) crystals located inside the coil of the magnet. These were used 
previously with great success by the CLEO II experiment. In BABAR and 
BELLE, another novel feature is a hadron calorimeter which is used for Kl 
and muon detection and allows the use of B^ ipKi in addition to the CP 
eigenstate B ip Kg. 
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Table 2: B Factory Detectors 


Name 

BaBar 

BELLE 

Logo 

Elephant 

Jazzy B 

^crossing 

0 mr 

±11 mr 

Vertex 

5 layers 

3 layers 

r (cm) 

3.2 - 14.4 

3-5.8 

CDC 

40 layers 

50 layers 

r (cm) 

24 - 80 

9-86 

PID 

DIRC+dE/dx 

Aerogel /TOE±dE/d; 

EM Cal 

CsI(Tl) 

CsI(Tl) 

Magnet 

1.5 T 

1.5 T 

Kp/fi 

RPC 

RPC 


Linseed Oil 

Glass 

DAQ 

Digital 

Q-to-T into 


Pipeline 

multihit TDC 


1.4 The First Measurements of CP Violation 

The measurement requires the reconstruction of fcp decays, the de¬ 

termination of the 6-flavor of the accompanying (tagging) B meson, the mea¬ 
surement of At, and a £t of the expected At distribution to the measured 
distribution using a likelihood method. 

BELLE reconstructs B^ decays to the following CP eigenstates: J/ipKs, 
ij{2S)Ks, XciKs, VcKs for ^ = -1 and for ^ = +!• BABAR 

includes all these modes except for rjcKg. The two classes = ±1) should 
have CP asymmetries which are opposite in sign. 

Both experiments also use B^ —>■ decays where K*^ 

Here the final state is a mixture of even and odd CP. The CP content can, 
however, be determined from an angular analysis of other ipK* decays. The 
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CP odd fraction is found to be small (i.e. 19 ±4 % (16 ± 3.5) in the BELLE 
(BABAR) analysis). 

The CP eigenstate event samples are shown in Figs, ^and The distri¬ 
butions of beam constrained mass are shown for the modes other than 
For the latter, the distribution of p*^, the momentum of the B candidate 
in the CM frame is shown. In the fully reconstructed modes, BELLE Ends 
747 events over a background of 58.6. For the 'ipKi sample, there are 569 
over a background of 223. The event yields from BABAR are similar. They 
observe a total of 1230 events over background of 200. It is clear that the 
CP eigenstate samples that are used for the CP violation measurements are 
large and clean. 

To identify the flavor of the accompanying B meson, leptons, kaons, 
charged slow pions from D* D^7r~^ decays, and energetic pions from two- 
body B decay (e.g. Zl*+7r+) are used. A likelihood based method 

is used by BELLE. BABAR first separately tags leptons and then kaons. 
If these methods fail or conflict, two neural nets are then used to iden¬ 
tify the remaining tags. The hgure of merit for flavor tagging is the ef¬ 
fective efficiency, Ce//, which is e(l — 2w‘^) summed over all tagging cate¬ 
gories. belle’s method gives Ce// = 0.270 ±0.008lo;oo9 while BABAR Ends 
Seff = 0.261 ± 0.012 for their method. 

The probabilities of an incorrect Eavor assignment are determined directly 
from the data using exclusively reconstructed, self-tagged ^ 

£){*)- 7 r+, D *~and decays. The values of wi are ob¬ 

tained from the amplitudes of the time-dependent B^B^ mixing oscillations: 
(A^^of — Nsf)/{Nqf -|- A"sf) = (1 — ‘^wi) cos(AmrfAf). Here Nqf and Nsf are 
the numbers of opposite and same Eavor events. These wrong-tagging rates 
determined from data are then used in the CP Et for BELLE. To incorporate 
this information in their analysis, BABAR performs a simultaneous Et to the 
Eavor tagged and CP eigenstate data samples. 

The fcp vertex is determined using lepton tracks from J/^jj or 'ip{2S) 
decays, or prompt tracks from pc decays. The /tag vertex is determined from 
well reconstructed tracks not assigned to fcp- Tracks that form a Ks are not 
used. For BELLE, the typical vertex-Ending efficiency and vertex resolution 
(rms) for zcp (^^tag) are 92 (91)% and 75 (140) pm, respectively. Note that 
the resolution on the tag side includes a large additional contribution from 
charm decay. For BABAR, the rms resolution in Az = zcp — Ztag is 180/rm 
and their vertexing efficiency is 97%. Once the tag and CP eigenstate vertices 
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are reconstructed, the proper time is calculated from 

Examination of the time dependent distributions from the two experi¬ 
ments (Figs. a 1.10 shows clear indications of CP violation. The raw time 
distributions for the and samples are indeed asymmetric. This is 
most striking in the BELLE hgure. Thus, CP violation is visible even in the 
raw time distributions. The time dependent asymmetry distributions for the 
= ±1 samples are indeed opposite as expected. Finally, the flavor tagged 
samples show no time dependent asymmetry. 

To extract the values of the CP violating parameter, the two experiments 
perform unbinned maximum likelihood hts to the time distributions of the 
tagged and vertexed events. The fits take into account the effects of back¬ 
ground, vertex resolutions and incorrect tagging. 

In the summer of 2001, the first signihcant measurements of the CP vi¬ 
olating parameter sin20i were obtained by BELLE and BABAR. BELLE 
hnds 

sin 201 = 0.99 ±0.14 ±0.06 

with a statistical significance of greater than six standard deviations^]. BABAR 
hnds 

sin 201 = 0.59 ±0.14 ±0.05 

with a signihcance of 4.1 standard deviations]^]. The results are based on 
data samples of comparable size (31 million and 32 million BB pairs, re¬ 
spectively). As discussed above, the efficiencies and resolutions of the two 
experiments are also quite comparable. However, the two measurements are 
only marginally consistent. Larger data samples and additional more precise 
measurements will be required to fully reconcile these two results. 


2 Determination of the other angles 

Experimental work on the determination of the other angles 02 and 03 has 
just started. To obtain precise results that test the validity of the Standard 
Model, much larger data samples will be required. Below I briehy describe 
some of the possible approaches. 

To measure 02 (a.k.a a), the two most promising approaches involve 
the use of the decay modes B^ —>■ and B^ The former 

is an example of CP eigenstate and is the most straightforward approach. 


The interference in this mode between the direct decay and the decay via 
mixing leads to a CP violating asymmetry as in the charmonium decay mode 

However, there are several additional complications. The decay ampli¬ 
tude for TT+TT” contains a contribution from a tree diagram (6 ^ uud) 

as well as a Cabibbo suppressed penguin diagram (6 —> sus). The penguin 
contribution is not negligible and has a weak phase that is different from the 
phase of the larger tree amplitude, which is zero in the usual parameteriza¬ 
tion. Therefore the time dependent asymmetry, proportional to sin(AmAt), 
which is measured is not equal to sin 202 but instead will have a large un¬ 
known correction. The presence of the extra contribution also induces an 
additional time dependent term proportional to cos(AmAf). This is called 
penguin pollution. 

There are a number of other purely experimental complications. The 
branching fraction for the decay is somewhat small. Using 

22.6 X 10® BB pairs BABAR hnds ^ tt+tt-) = (4.1 ± 1.0 ± 0.7) x 10“® 
[y], while using a sample of 11.1 x 10® BB pairs BELLE hnds B{B^ 
7r~^7T~) = (5.6 ± 2.3 ± 0.4) x 10“® |^. Moreover, even after the application 


of high momentum particle identihcation, the CP eigenstate signal sits on a 
large continuum background. The signal to continuum background ratio is 
a function of the tag method; this effect must be taken into account in the 
CP extraction. The conventional continuum suppression variables are also 
correlated with tagging. There is still residual background from misidentihed 
R® ^ K^tt~ as well. 

Despite these difficulties, the hrst attempts to determine sin 202 using 
the 7r’''7r“ mode have started. BABAR reports results of a preliminary £t 
which gives sin202(e//) = O.O30o;fg ± 0.11. BELLE expects a sensitivity of 
±0.6 and will present a result soon. The size of the correction for penguin 
pollution in sin 202 (e//) is unknown. 

The angle 03 is particularly difficult to measure. There are at least three 
methods which give some information on this angle. 

The method which requires the least data but has the largest theoretical 
uncertainty involves the measurement of the branching fractions of the Ktt 
and TTTT decay modes. The underlying idea is to compare decays where the 
interference of the penguin and tree amplitudes are different to extract about 
information about the phase 03. For example, if 03 > 90®, then the interfer¬ 
ence term is constructive for R® — K^tt~ and destructive for R® —> tt+tt". 
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leading to a ratio for the widths of the two modes greater than one. Con¬ 
versely, if 03 < 90*^ then a K~'k ^ratio less than unity is expected. 
The current data somewhat favor 03 > 90° whereas the indirect CKM hts 
prefer 03 ~ 60°. The interpretation of this discrepancy is still unclear. 

A theoretical industry has been developed to extract the maximum pos¬ 
sible information on 03 from these and related decay modes|^. However, it 
has been argued by some theorists that other smaller diagrams or hnal state 
interaction effects may invalidate these analyses [p!^. 

Another method to determine 03 uses the Cabibbo suppressed decay mode 
B~ —> In the case that fcp, where fcp is a CP eigenstate 

such as K~K~^ , KsTt^, two amplitudes with different weak phases will 

interfere. The asymmetry and rate again involve the phase of Vub- This 
method is theoretically cleaner than B Ktt but requires extremely large 
data samples. The current yield of B~ —> DcpK~ events from BELLE is 
~ 16 events in a 21 fb“^ data sample |]T 6 | . This yield may be increased by 
about a factor of two by including additional CP eigenstates of the D meson. 
Another approach is to use doubly Cabibbo suppressed decays (DCSD) of 
the D meson. In this case, the event yield is much smaller but the direct CP 
asymmetries are much larger^ISf. 

Measurement of a time dependent asymmetry in the mode 5° ^ 
also gives information on the angle 03. In this case, however, the asymmetry 
is proportional to sin(20i -|- 03). The time dependence is also more compli¬ 
cated than in the case of a simple CP eigenstate such as i?° ^ To 

reconstruct 5° ^ , one uses a partial reconstruction technique where 

the 7 r“ and the low momentum pion from the D*^ decay are 

observed but the Zl° meson is not reconstructed. This increases the number 
of useable signal candidates by roughly an order of magnitude. Initial MC 
estimates indicate that, however, that even with the partial reconstruction 
method, 200 fb“^ (which corresponds to 17000 lepton tagged signal events) 
will give a precision of ±0.34 for sin(20i ± 03) in this channel [P!^ . 

For the methods described above one usually measures a trigonometric 
function of one of the angles such as sin(20i). In such a case, there often 
exists multiple solutions or ambiguities for the angle itself (e.g. four solutions 
for 01 from the initial sin(20i) measurements). One cannot simply discard 
the solutions that are not consistent with the Standard Model. Instead to 
resolve these ambiguities, requires difficult high statistics measurements of 
quantities that depend on other functions of the angles such as cos( 0 ). 
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In addition to the program of measnring the angles of the nnitarity trian¬ 
gle, there is also the qnestion of whether there are new CP violating phases 
from new interactions or physics beyond the Standard Model. At the mo¬ 
ment, such new phases are completely unconstrained. One way to attack 
this question is to measure the time dependent CP asymmetry in penguin 
modes such as (pKg or ^ rj Kg and compare it to the asymmetry in 

B^ 'ipKg. In the absence of new physics, they should be equal. However, 
if there are new physics contributions in penguin loops these asymmetries 
will differ substantially. Again an order of magnitude more data is needed 
for stringent tests|^,[^. This search for new physics in CPV will be one 
of the most interesting aspects of the next phase in B factory physics. 
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Figure 1: BABAR CP eigenstates. The upper figure shows the beam con¬ 
strained mass distribution of the fully reconstrructed CP eigenstates. The 
lower figure shows the distribution for the —> 'ijjKi candidates. 
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Figure 2: Belle CP eigenstates: the upper figure shows the beam constrained 
mass distributions for the exclusive modes while the lower figure shows the 
p% distribution for the 'iPKl candidates. 


14 
























-5 0 5 

At (ps) 


Figure 3: BABAR data on time dependent CP asymmetry. For the fully 
reconstructed CP = — 1 final states, (a) The At distributions for tags; 
(b) The At distributions for tags; the shaded area is the background, the 
solid curve is the £t to the data, (c) The time dependent asymmetry between 
(a) and (b); For the 'ipK^ final state, (d) The At distributions for B^ tags; 
(e) The At distributions for B^ tags; (f) The time dependent asymmetry 
between (d) and (e); 
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Figure 4: BELLE data on time dependent CP asymmetry. At distributions 
for the events with q^f = +1 (solid points) and q^f = —1 (open points). The 
results of the global £t (with sin2(/)i = 0.99) are shown as solid and dashed 
curves, respectively. 
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Figure 5: BELLE data on time dependent CP asymmetry. The raw A dis¬ 
tributions for and tags, (a) The asymmetry obtained from separate 
hts to each At bin for the full data sample; the curve is the result of the 
global £t. The corresponding plots for the (b) {cc)Ks (^/ = —1), (c) 

(^/ = and (d) B^ control samples are also shown. The curves are the 
results of the £t applied separately to the individual data samples. 
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